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The photoionization cross section a, spin-polarization pa- 
rameters P and Q, and the angular-distribution asymmetry 
parameter /3 are calculated for the 7s state of francium for 
photon energies below 10 eV. Two distinct calculations are 
presented, one based on many-body perturbation theory and 
another based on the model potential method. Although pre- 
dictions of the two calculations are similar, the detailed en- 
ergy dependence of the photoionization parameters from the 
two calculations differ. From the theoretical p-wave phase 
shifts, we infer quantum defects for p 1 / 2 and P3/2 Rydberg 
series, permitting us to calculate positions of experimentally 
unknown p states in francium. 
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I. INTRODUCTION 

Remarkable progress has been made recently in de- 
termining energies and lifetimes of low-lying states of 
the heaviest alkali-metal atom francium jjj, motivated 
in part by the enhancement of parity non-conserving 
(PNC) effects in francium compared with other alkali- 
metal atoms. This experimental work has been accom- 
panied by theoretical studies of properties of the francium 
atom [QgJ , concerned mostly with energies and hyperfine 
constants of the ground and low-lying excited states or 
transitions between such states. 

In this work, we present two calculations of photoion- 
ization of francium for photon energies below 10 eV; the 
first is an ab-initio many-body calculation and the sec- 
ond is a model potential (MP) calculation. Experiments 
on photoionization of francium are planned for the Ad- 
vanced Light Source at Berkeley. Q 

Ab-initio calculations of photoionization in alkali-metal 
atoms have proved to be a formidable challenge. Pho- 
toionization calculations in cesium based on the Dirac 
or Breit-Pauli equations accounted for the spin- 

orbit interaction, but not for shielding of the dipole op- 
erator by the core electrons or for core-polarization ef- 
fects; whereas, relativistic calculations that included cor- 
rections from many-body perturbation theory (MBPT) 
at the level of the random-phase-approximation (RPA) 
accounted for the spin-orbit effects and for core 



shielding, but not core polarization. Predictions from 
these many-body calculations were in poor agreement 
with measurements of the Fano spin-polarization param- 
eter P by Heinzmann et al. [|10|1, of the spin-polarization 
parameter Q by Lubell et al. [[1 if , and with the measure- 
ment of the angular-distribution asymmetry parameter (3 
by Yin and Elliott The first quantitatively success- 
ful many-body calculation of photoionization of cesium 
was a relativistic many-body calculation that included 
both core polarization and core shielding corrections [^3) ; 
that method is applied to low-energy photoionization of 
francium in the present paper. 

Although successful many-body calculations of pho- 
toionization of heavy alkali- metal atoms are of recent vin- 
tage, nearly three decades ago, a number of increasingly 
sophisticated and successful model potential calculations 
of the photoionization of cesium appeared |l4|-p"7|] , cul- 
minating with that of Norcross p8[ . The latter calcu- 
lation, which included the spin-orbit interaction, long- 
range polarization potentials, and shielding corrections 
to the dipole operator, gave quantitatively correct val- 
ues for all of the measured photoionization parameters 
in cesium. A model potential similar to the one used 
in 1 18 1 was developed recently |3) to study transitions in 
francium and is used here to study low-energy photoion- 
ization in francium. 

Below, we sketch the important features of the theo- 
retical methods. The photoionization cross sections for 
Fr 7s are calculated in Section ILB in both methods and 
are compared against one another. In Section II. E, we 
give results for the spin-polarization parameters and the 
angular distribution asymmetry parameter. Section III 
concludes our discussion of the francium photoionization. 



II. THEORETICAL ANALYSIS & DISCUSSION 

A. Many-Body Perturbation Theory 

We start our many-body analysis from the Dirac- 
Hartree-Fock (DHF) Vn-i approximation, in which the 
DHF equations are solved self-consistently for core or- 
bitals, and the valence orbitals are determined subse- 
quently in the field of the "frozen core" . The total phase 
shift 5 K for a continuum state with angular quantum 
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number k in the field of the core is a sum of rapidly vary- 
ing Coulomb phase shift <5j? and the short-range shift 6 K . 
The short-range DHF phase-shifts for pyn and p 3 / 2 con- 
tinuum wave functions are shown in Fig. The p 3 / 2 wave 
function lags in phase compared to the P1/2 wave func- 
tion owing to the spin-orbit interaction, which is attrac- 
tive for pi/2 states and repulsive for p 3 / 2 states. The DHF 
approximation typically underestimates removal energies 
of bound electrons in heavy atoms such as francium by 
about 10%; similar accuracy is expected for phase shifts. 
To improve this level of accuracy one must take into ac- 
count higher-order MBPT corrections. 

The clear advantage of the Vjq-i approximation stems 
from the fact that one-body contributions to the residual 
Coulomb interaction vanish. This leads to a significant 
reduction in the number of terms in the order-by-order 
MBPT expansion. In particular, first-order corrections 
to the energy (or the phase shift) vanish and the pertur- 
bation expansion starts in second-order. 

The leading correlation contribution to the energy is 
the expectation value of the second-order self-energy op- 
erator given diagrammatically by the Brueckner- 
Goldstone diagrams of Fig. |[ Solutions to the Dirac 
equation including the Vn-i potential and the self-energy 
operator are called Brueckner orbitals (BO). The non- 
local self-energy operator S, in the limit of large r, de- 
scribes the interaction of an electron with the induced 
electric moments of the core, 



Z tj (r) = 1 + (z - l)e- a i' r + r(a| i) + 



S(r,r',e) 



2r 4 



5{r-r') 



(1) 



where is the dipole polarizability of the core. We 
determine second-order correction to the phase shift per- 
turbatively as 



5^ =-sin- 1 (7r (u eK \^\u eK )) 



Here, u eK is a continuum DHF wave function normalized 
on the energy scale. The resulting DHF+BO phase shifts 
are presented in Fig. |l|. The attractive polarization po- 
tential draws in the nodes of the wave function, resulting 
in larger phase shifts. The change in the phase shift is ap- 
proximately the same for both p 3 / 2 and P1/2 continuum 
states, demonstrating that the self-energy correction is 
mainly due to the accumulation of phase outside of the 
core. 



B. Model Potential 

The parametric model potential used in this work has 
the form pi) 



r 2r 4 1 



= _( r/r O))6 



(3) 



where ad is the static dipole polarizability of the Fr + 
ionic core and the effective radial charge Zij(r) is given 
by 



4 J V)e- Q 2 



(4) 



The angular momentum-dependent parameters, a f , i 

1,...,4 and the cut-off radius Tc are obtained through 
a non-linear fit to one-electron Rydberg energy levels in 
francium |l],^,[l9). Because the spin-orbit effects are ap- 
preciable for heavy alkali metals, two separate nonlinear 
fits; one for each fine-structure series, j+ = I + | and 
j_ = I — 5 were performed. The static dipole polarizabil- 
ity was obtained from an extrapolation of the known core 
polarizabilities for the other alkali metals as a<j(0) = 23.2 
a.u. ||]. We note that an ab initio value for the francium 
core polarizability is now available |2j . 

A comparison of short-range phase shifts calculated in 
the model-potential method and the MBPT is presented 
in Fig. ||. We find reasonable agreement between the two 
methods. The MP continuum wavefunctions are slightly 
lagging in phase compared to many-body wavefunctions. 
Such phase differences result in Cooper minima being 
shifted to higher photoelectron momentum in the model- 
potential calculation. 



C. Quantum defects 

In quantum defect (QD) theory ^(|, the energy levels 
of the valence electron are described by a hydrogen-like 
Rydberg-Ritz formula 
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2(n-/i K ) 2 



(5) 



in terms of a quantum defect p Kl which is represented 
as an expansion in powers of energy with constant coef- 



(2) ficients p} K 



u -u {0) +u (1) e +u [2) (e ) 2 + 



(6) 



The Rydberg-Ritz formula provides an accurate fitting 
expression for the bound spectrum of alkalis. The QD 
is related to threshold value of the phase shift as 
fj,n = 5 k /tt + integer. The QD's for Fr p states are not 
known, since the relevant Rydberg series have not been 
observed experimentally. We use our ab-initio threshold 
phase shifts together with experimentally known energies 
for 7p and 8p states to predict QD's; thereby approxi- 
mating the entire Rydberg spectrum of Fr p states. The 
predicted quantum defects are given in Table |. We as- 
signed an error bar of 0.5% to the threshold phase shift, 
based on the accuracy of an application of the many-body 
formalism employed here to the case of Cs Q . 

In Table |[ we also present the MP values of quantum- 
defects obtained by fitting Rydberg series calculated with 
the potential in Eq. E| We find generally good agreement 
for the leading order quantum-defect /i« , estimated in 
the two methods. Higher-order QD parameters, /i^ and 

(2) 

p K , calculated in the two methods do not agree well. 
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This is due to the sensitivity of these parameters to the 
value of //« . The values for /i« obtained by fitting to 
the MP-calculated np levels agree to four significant dig- 
its with the values, for /il°' > , extracted from the threshold 
phase shifts in Fig. ||. 

Using the calculated quantum defects, we predict en- 
ergy levels for the lowest few np states. Table || lists 
these energies and compares them with the present MP 
calculation and with a recent MBPT single-double (SD) 
calculation j22| . The accuracy of our many-body calcula- 
tion was estimated by exercising upper and lower bounds 

on and a consistent determination of fJlp and 
to fit 7p and 8p energies. MBPT results are in reasonable 
agreement with the MP calculations and SD predictions 
for these levels. 



D. Cross-section 

The total cross-section for photoionization of the va- 
lence electron v is the sum of partial cross-sections 

K K 

where u> is the photon energy. The dipole transition am- 
plitude for an ionization channel v — > en is defined as 

D« = r m e^<K eK ||r|K), (8) 

where u v is the valence wave function and where the con- 
tinuum wave function u tK is normalized on the energy 
scale. Here we have two ionization channels 7s — > ep±/2, 
with k = 1 and 7s — > £^3/2, with k — —2. The DHF 
results for the total cross-section are shown with dashed 
lines in Fig. ^. Since the DHF potential is non-local, the 
resulting amplitudes depend on the gauge of the elec- 
tromagnetic field. The difference between length- and 
velocity-form values is especially noticeable in the near- 
threshold region. 

Second-order corrections, and the associated all-order 
sequence of random-phase approximation (RPA) dia- 
grams, account for the shielding of the external field by 
the core electrons. Explicit expressions for the second- 
order MBPT corrections can be found, for example, in 
Ref. Already in second order, the dipole operator 

with RPA corrections reduces at large r to an effective 
one-particle operator 

where ay (<•<-> ) is a dynamic polarizability of the core. The 
first term is associated with the applied electric field and 
the second with the field of the induced dipole moment of 
the atomic core; the valence electron responds to a sum 
of these two fields. We note that the induced field may 
become strong and reverse the direction of the total field. 



The RPA cross-section is presented with a thin solid 
line in Fig. ||. In contrast to DHF amplitudes, the 
RPA amplitudes are gauge- independent. Furthermore, 
we note the sudden upturn in the RPA cross-section for 
the photoelectron momenta p rs 0.7 a.u. associated with 
a J = 1 core excitation resonance. To predict the posi- 
tion of this resonance, we calculate the dynamic polariz- 
ability of Fr + within the framework of relativistic RPA, 
discussed in p4| ). The energy of the first resonance is at 
Lu r = 0.4024 a.u.. Using the DHF value of 7s threshold, 
0.1311 a.u., we expect the first core excitation resonance 
to appear at p rs 0.74 a.u.. The dynamic polarizability 
of Fr core ced(u>) from this RPA calculation is plotted as 
a function of electron momentum p in Fig. ^. 

To account for core-polarization corrections to the 
DHF wave function, discussed in the introduction, we 
evaluate the second-order corrections to the DHF wave 
functions of the valence electron due to the self-energy 
operator E^ 2 ^ 

V (2) 
(2) 

The resulting orbital u v +u v is the perturbative approxi- 
mation to the valence-state Brueckner orbitals (BO). Ap- 
proximate Brueckner orbitals for a continuum state (e k) 
are found by solving the inhomogeneous Dirac equation 

(h + Vtf-i-e) w eK = (-7rsin(5 K -I] (2) ) u tK (11) 

normalized on the energy scale, where S K is given in 
Eq. (||). Brueckner orbitals for the 7s valence state and 
a P1/2 continuum state are compared with unperturbed 
DHF orbitals in Fig. |. 

The BO corrections contribute to transition ampli- 
tudes starting from third order. Together with the RPA 
corrections, they provide the most important third-order 
contributions for bound-bound transitions, as discussed 
in [ p5| . In the present approach, we modify the conven- 
tional RPA scheme by replacing the valence and contin- 
uum wave functions by the approximate Brueckner or- 
bitals described above (RPA0BO). Such a modification 
accounts for the important second- and third-order cor- 
relation corrections and for a subset of fourth-order con- 
tributions to transition amplitudes. We note that this 
fourth-order subset brings the photoionization parame- 
ters in cesium into good agreement with available exper- 
imental data; therefore, we believe that this approach 
will provide reliable predictions for francium. The re- 
sulting cross section is shown with a heavy solid line 
in Fig. ^, and decomposed into partial cross-sections in 
Fig. 0. Calculations using length and velocity forms of 
transition operator lead to slightly different result in the 
modified RPA0BO scheme; we present the final result 
in the length form only. Both photoionization channels 
exhibit Cooper minima; o~ Pl/2 vanishes at p w 0.1 a.u. 
and (Jp 3/2 vanishes at p w 0.5 a.u.. Combining the two 
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partial cross sections, leads to a broad minimum in the 
total cross-section slightly below p = 0.45 a.u.. The total 
cross section in Fig. is not very sensitive to the posi- 
tions of Cooper minima in the P1/2 and P3/2 channels. 
Conversely, the spin-polarization and angular distribu- 
tion measurements, discussed in the following section, 
provide information sensitive to fine details of individual 
transition amplitudes. 

Fig. H examines the total photoionization cross sections 
for francium, calculated in the two method. The label 
" static" refers to the set of MP results with the core static 
dipole polarizability in Eq. ||. The shielding of the elec- 
tron dipole operator is truncated in the MP calculations 
by introducing a cut-off term, similar to the exponential 
term in the one-electron potential in Eq. The threshold 
cross sections in the P1/2 and P3/2 channels (not shown 
here) are, respectively, 1.74 and 0.02 Mb. Cooper minima 
appear in both channels at approximately, p w 0.15 a.u. 
and p « 0.75 a.u. and the maximum cross section in the 
P1/2 channel is cr Pl/2 (max) w 0.2 Mb. The Cooper min- 
imum in the P3/2 photoelectron cross section calculated 
in the MP method with the static core dipole polarizabil- 
ity, occurs approximately where the first core resonance 
in Fig. becomes excited. By including the dynamic core 
polarizability otd{oj) in the MP calculations, the curve la- 
beled as "dynamic" in Fig. || is obtained. The Cooper 
minima are moved to lower photoelectron momenta, re- 
sulting in a shallow minimum in the total cross section 
near p w 0.5 a.u. 

The comparison in Fig. |§| indicates that the cross sec- 
tions calculated in the MP method are in general larger 
than the MBPT cross sections. The "MP dynamic" and 
the MBPT cross sections both rise for the photoelectron 
momenta p > 0.5 a.u. to meet the first core-excited res- 
onance near p ~ 0.75 a.u. 



E. Polarization parameters 

Fano jl0| proposed a measurement of spin polarization 
P of photoelectrons emitted from unpolarized Cs atoms 
illuminated by circularly polarized photons. The total 
spin polarization is expressed in terms of P1/2 and P3/2 
transition amplitudes as (23| Q 

5 \D 3/2 \ 2 - 2 \D 1/2 \ 2 + 4^2 X[D 1/2 D* /2 ] 

6(\D 3/ 2\ 2 + \D 1/2 \ 2 ) ' (1 } 

The result of our RPA0BO calculation of the spin- 
polarization parameter P is presented in Fig. 0, where 



it is seen that the polarization reaches 100% at momen- 
tum p w 0.3 a.u. The model-potential results for P are 
also given in Fig. |9[and compare well with the RPA©BO 
calculation in FigT^j. Maximum spin polarization in the 
MP method occurs at p ~ 0.35 a.u. The calculations 
with the static and dynamic core polarizabilities in Eq. ^ 
are similar and differ only after the maximum is reached. 

Lubell and Raith [[LlJ measured a different spin- 
polarization parameter Q obtained from photoionization 
of polarized Cs atoms by a circularly polarized light. In 
the Lubell- Raith setup, the p 3 / 2 channel can be accessed 
individually, for example, by photoionization with lcft- 
circularly polarized light of the 7s, electron prepared in 
the m s = +5 substate. Combining the partial cross sec- 
tion <7 P3/2 thereby obtained with the total cross-section, 
permits one to deduce the partial cross-section for the 
P1/2 channel. The Lubell- Raith parameter Q is defined 
as the ratio of the difference to the total of the photoab- 
sorption intensities for two photon helicities 



Q 



\D3/2\ 2 -2\D 1/2 \ 2 
2(\D 3/2 \^ + \D 1/ 2\ 2 ) 



(13) 



The limiting values for the Lubell-Raith parameter are 
— 1 < Q < h- We stress that a measurement of Q or 
of the (phase-insensitive) parameter P, together with a 
measurement of the total cross-section permits one to ob- 
tain information about absolute values of transition am- 
plitudes. A further measurement of the phase-sensitive 
angular-distribution parameter (3 |23], 



\D 3 ,2? 



2V25R 



D 1/2 D 



3/2 



\D 3/ 2\ 2 



\D lf 2\ 2 



(14) 



would permit one to determine the relative phase be- 
tween the P3/2 and pi/ 2 continuum amplitudes and would 
constitute an essentially complete description of the pho- 
toionization process. The many-body result for (3 is 
shown in Fig. ^. The differential cross section is pro- 
portional to 1 — i/3P2(cos#). 

The MP results for Q and the asymmetry parame- 
ter (3 are given also in Fig. [5| The comparison between 
MP and MBPT results is generally favorable: the results 
with the core dynamics polarizability in Eq. |9| are in bet- 
ter qualitative agreement with the MBPT calculations. 
We note that near p w 0.45 a.u., the photoelectron has 
the propensity to be ionized perpendicular to the photon 
polarization axis and near p « 0.6 a.u., the photoelec- 
tron is preferentially ejected in the j ' = k channel, where 
Q — > — 1. A similar situation is evident from MBPT re- 
sults at p ~ 0.5. The other limiting value is reached near 



threshold, where a 



Pl/2 



There is a phase difference in the _D 1/ / 2 £ ) 3/2 interfer- 
ence term in Eqs. [j^,[t4] and the corresponding equations in 
Ref. caused by the unconventional definition of reduced 
matrix elements used in that work. 



III. CONCLUSION 

We have calculated the photoionization cross sections 
of the ground-state francium. Both many-body and 
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model-potential approaches were employed to obtain the 
cross sections, quantum defects, spin-polarization param- 
eters and photoelectron asymmetry parameter. We find 
Cooper minima in both^/2 andp 3 / 2 channels. The com- 
parison between the MBPT and MP results are satisfac- 
tory. The Cooper minima predicted in the MP calcu- 
lations are at higher photoelectron energies than those 
calculated in the MBPT method. The origin of this 
difference can be traced to the shielding of the valence- 
electron dipole by the core electrons. The induced dipolc 
moment of the core manifests itself as a dynamic polar- 
izability term. Upon replacing the static core polariz- 
ability with the dynamic polarizability, better quantita- 
tive agreement with the MBPT results is observed. We 
predict the energy-dependence of the photoelectron spin- 
polarization and asymmetry parameters which we hope 
will stimulate further experimental work in francium. 
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FIG. 4. Total photoionization cross-sections for Fr 7s 
state, calculated in various many-body approximations. 
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FIG. 1. Short-range phase shifts for p continuum in Fr cal- 
culated in Dirac-Hartree-Fock (DHF) approximation and in- 
cluding self-energy correction (BO). 
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FIG. 5. RRPA dynamic polarizability of Fr + as a func- 
tion of photoelectron momentum, calculated with the 7s DHF 
threshold, 0.13107 a.u.. 
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FIG. 2. The second-order self-energy operator. 
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FIG. 3. Comparison of short-range phase shifts obtained in 
model potential and many-body methods. Upper two curves 
represent pi/2-phase shifts, and two lower curves represent 
those for p 3 / 2 continuum. 
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FIG. 6. Brueckner (BO) and DHF orbitals. Upper panel: 
large component of continuum wave function; Lower panel: 
radial density. 
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FIG. 7. Partial and total cross-sections for Fr 7s calcu- 
lated in the RPAffiBO many-body approach. 
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FIG. 8. Comparison of Fr 7s photoionization cross sec- 
tions calculated with the MBPT and the MP methods. The 
curves labeled as "static" and "dynamic" refer, respectively, 
to the calculations with static and dynamic core polarizabili- 
ties in Eq. []. 
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FIG. 9. Fano P and Lubell-Raith Q spin-polarization pa- 
rameters, and dipole asymmetry parameter (3, calculated in 
the RPA0BO (solid line) and MP approach. The long-dashed 
curve represents the MP results using the RRPA dynamic po- 
larizability in the transition operator, and dotted curve — 
static polarizability instead. 



TABLE I. Predicted quantum-defect parameters 



state 


u (0) u (1) 






MBPT 




Pl/2 


4.620(3) -0.3(1) 


5.3(7) 


PZ/2 


4.517(3) -0.5(1) 


4.7(8) 




Model Potential 




Pl/2 


4.605 -0.79 


1.61 


PZ/2 


4.505 -0.87 


1.59 
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TABLE II. Predicted p levels of Fr in cm -1 . 



Present Work 





odt+mrpt 3 


MP 


ST) b 


^i J l/ 2 


5748(3) 


5737 


5738 


10pi/2 


3800(2) 


3790 


3795 


HPl/2 


2700(2) 


2691 




12pi/ 2 


2016(1) 


2011 




9P3/2 


5496(2) 


5487 


5488 


10p 3 /2 


3662(2) 


3655 


3659 


11P3/2 


2616(2) 


2610 




12p 3 /2 


1962(1) 


1958 





a The error bars were estimated based on the agreement [[13| of 
the predicted second-order phase shifts with known quantum 
defects for cesium p-states. 
b Values marked as predicted in Ref. p2j. 
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